Dynamics of a vortex domain wall in a magnetic nanostrip: an application of the 

collective-coordinate approach 
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The motion of a vortex domain wall in a ferromagnetic strip of submicron width under the influence 
of an external magnetic field exhibits three distinct dynamical regimes. In a viscous regime at low 
fields the wall moves rigidly with a velocity proportional to the field. Above a critical field the viscous 
motion breaks down giving way to oscillations accompanied by a slow drift of the wall. At still higher 
fields the drift velocity starts rising with the field again but with a much lower mobility dv/dH 
than in the viscous regime. To describe the dynamics of the wall we use the method of collective 
coordinates that focuses on soft modes of the system. By retaining two soft modes, parametrized 
by the coordinates of the vortex core, we obtain a simple description of the wall dynamics at low 
and intermediate applied fields that describes both the viscous and oscillatory regimes below and 
above the breakdown. The calculated dynamics agrees well with micromagnetic simulations at low 
and intermediate values of the driving field. In higher fields, additional modes become soft and 
the two-mode approximation is no longer sufficient. We explain some of the significant features of 
vortex domain wall motion in high fields through the inclusion of additional modes associated with 
the half-antivortices on the strip edge. 
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I. INTRODUCTION 



Dynamics of domain walls in ferromagnetic strips and 
rings with submicron dimensions is a subject of active 
research J^i^ This topic is directly relevant to several pro- 
posed schemes of magnetic memory and is also interesting 
from the standpoint of basic physics. The dynamics of 
domain walls under an applied magnetic field has distinct 
regimes: viscous motion with a relatively high mobility 
at low fields and underdamped oscillations with a slow 
drift at higher fields 4 

In a nanostrip, domain wall dynamics is further com- 
plicated by the composite nature of the wall which con- 
sists of a few ~ typically two or three ~ elementary topo- 
logical defects in the bulk and at the edge of the strip.- As 
a result, its motion is dominated by a few low-energy de- 
grees of freedom associated with the motion of the topo- 
logical defects. Weak external perturbations engage only 
the softest (zero) mode - a rigid translation of the do- 
main wall along the strip. Larger external forces excite 
additional modes thereby altering the character of mo- 
tion. 

The general approach to the dynamics of composite 
domain walls was described recently by Tretiakov et al^ 
The configuration of a domain wall is parametrized by a 
few collective coordinates ^ — {^i, ^2, • ■ • , ^n} represent- 
ing the soft modes, and the free energy of the system U is 
treated as a function of ^. The Landau-Lifshitz equation 
for the spin dynamics with damping in Gilbert's form^'* 
is translated into a set of coupled equations of motion 
for the collective coordinates. In a vector notation, they 



-1/2 



i ■ 



\ i W >^ '•■^ N S N N- 

.^. \ \ \ V, ^ \^ . ,' , \ \ V- 



-1/2 



FIG. 1: (Color online) A head-to-head vortex domain wall in 
a long strip of permalloy 200 nm wide. The thickness in the 
out-of-plane direction is 20nm. Simulation using OOMMF.— 
The winding numbers of the three topological defects are la- 
beled. The shaded region indicates out-of-plane magnetiza- 
tion within the vortex core. 



read 



F - + G| = 0. 



(1) 



Here components of the vector F are generalized forces 
Fi — —dU/d£,i derived from the free energy U; the 
symmetric matrix Fy characterizes viscous friction; the 
antisymmetric gyrotropic matrix Gij describes a non- 
dissipative force of a kinematic origin related to spin pre- 
cession. 

In this paper we describe in detail the application of 
the collective-coordinate method to the motion of a par- 
ticular type of domain walls, expanding on our previous 
report.^ Specifically, we focus on vortex domain wallsi^ 
found in long permalloy strips of submicron width w and 
thickness t <^ w. In such a strip, the magnetization pat- 
tern is essentially uniform across the thickness, allowing 
us to consider patterns that vary only along the length 
(x-direction) and width (y-direction) of the strip. 



2 



A vortex domain wall is a magnetization pattern that 
consists of three topological defects (Fig.[T|): a vortex in 
the bulk of the strip and two half-antivortices confined 
to the edgesjii We will show that the seemingly complex 
dynamics of a domain wall can be reduced to a simple 
motion of these defects. The collective-coordinate ap- 
proach focuses on the soft modes of a system. In the 
case of a vortex domain wall, the two softest modes turn 
out to be the coordinate {X) of the vortex core along the 
long axis of the strip and the coordinate (Y) of the core 
across the strip width. The dynamic properties of the 
vortex are similar to that of a massless charged particle 
moving through a viscous medium in the presence of a 
potential U {X, Y) and a fictitious magnetic field directed 
normal to the plane. The "electric charge" of the parti- 
cle equals Att times the topological charge of the vortex,= 
while the strength of the "magnetic field" equals the two- 
dimensional spin density. 



At low fields the equations describe steady viscous mo- 
tion of the wall with a velocity proportional to the applied 
field. The vortex is shifted in the transverse direction by 
an amount proportional to the velocity of the wall. At a 
certain critical field the vortex is expelled from the strip 
and the steady motion breaks down, giving way to an 
oscillatory regime where the vortex periodically crosses 
the strip. Each time the vortex comes to an edge it is 
expelled and reinjected to cross in the opposite direction. 
At very large magnetic fields one reaches an extreme os- 
cillatory regime, where the vortex quickly moves back 
and forth across the strip almost along the equipotential 
lines, while slowly drifting in the direction of the applied 
field due to viscous forces. 



The comparison of the predicted velocity curve to 
experimental^ and numerical^ data shows that the the- 
ory agrees well with the observations for low fields and 
fields just above the breakdown. However, the theory 
predicts a significantly lower drift velocity at high fields 
than is observed. The origin of the discrepancy has been 
traced to the appearance of additional soft modes in this 
regime that are not taken into account in the two-mode 
approximation. The new modes increase the amount of 
dissipation in the wall and thereby lead to faster drift. 



The paper is organized as follows. In Section |TT] we re- 
view the general aspects of the collective-coordinate ap- 
proach in applications to the dynamics of magnetization. 
In Section IIIII we derive the equations of motion for the 
two softest modes of the vortex wall parametrized by the 
coordinates of the vortex core and discuss the general 
aspects of their dynamics. A detailed analysis of the dy- 
namics is given in Section [TVl The additional soft modes 
are discussed in Section |Vl Auxiliary results are derived 
in the Appendixes. 



II. COLLECTIVE COORDINATES 

A. General formalism 

The dynamics of magnetization M(r, t) in a ferromag- 
net well below the Curie temperature is described by the 
Landau-Lifshitz-Gilbert equation^'^ for the unit vector 
m = M/M, 

dm dm 

— = 7^oHoff X m -I- am X — , (2) 

where /ioHofr(r) = —6U/S'M.{r) is an effective magnetic 
field. The gyromagnetic ratio 7 — g\e\/{2me) is 1.75 x 
10^^ s Akg^-'^ and the Gilbert damping constant a ~ 0.01 
in permalloy.iS 

The free energy includes, at the very least, exchange 
and dipolar interactions as well as the Zeeman energy of 
the ferromagnet in an external field H, 



U = JdV (A| Vm|2 + - a^qH • M) , 



(3) 



where Hin is the field induced by the nonuniform mag- 
netization. It satisfies equations V • (Hi„ -I- M) — and 
V X Hin = 0. The exchange constant in permalloy is 
A ~ 1.3 X lO^^'^ J/m. Crystalline anisotropy is negligibly 
small in permalloy; however, such terms can be included 
in the free energy ([3|). 

In principle, an infinite number of coordinates are nec- 
essary to describe the time evolution of a magnetization 
texture. However, as a domain wall propagates along 
a nanostrip, much of its structure remains unchanged. 
While the wall may distort, it retains (at least temporar- 
ily) its general shape. For instance, in a vortex wall the 
chirality of the vortex does not change while the vortex 
remains within the strip. Likewise, the magnetization far 
from the wall remains fixed as the wall moves. This sug- 
gests that the motion of the wall may be described by 
a finite set of coordinates ^ = {Ci:'C2: • • ■ ,£,n}, so that 
m = m(^(i),r). In particular, the evolution of magneti- 
zation is related to changes in the generalized coordinates 
as 

dm 9m • , ^ 

-df - Wfi^ 

where a sum over the repeated indices is implicit. By 
substituting Eq. (|4]) into Eq. ([2]) and integrating over 
the volume of the sample, we arrive at the equations of 
motion^ 

- r,,4- + G,,ij - 0, (5) 

with antisymmetric gyrotropic matrix Gy and symmetric 
matrix of viscosity coefficients . Here 

5U dm _ dU 
8m d£,i ' 

9m dm 



F, 



dV 
aJ f dV 



. dm dm 

J I dV m- -— X -— , 

d^i d^j 



(6a) 
(6b) 
(6c) 
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and J = M/7 is the density of angular momentum. 
Note that Gy obeys the identity 

9Gij dGjk dGkt ^ , . 

If the space of collective coordinates is simply connected, 
then one may express the gyrotropic tensor in terms of 
a gauge field Af. dj = {dAj/d^i) - {dAJd^j). The 
equations of motion ^ may then be derived from the 
Lagrangian 

L - A,i, ~ {/, (8) 
together with the Rayleigh dissipation function R = 

B. Soft and hard modes 

Equations (O-® are formally exact when they take 
into account all of the modes of a magnetic texture. If we 
are not interested in such level of detail, we may focus 
on a few modes that capture the most salient features 
of magnetization dynamics and neglect all other modes. 
It is useful to divide modes into soft ones, which remain 
active on the typical time scale T of the dynamics, and 
hard ones, the motion of which decays on a much shorter 
time scale. Since the drift velocity of a domain wall in 
an applied magnetic field is ultimately determined by 
the rate at which its Zeeman energy is dissipated, soft 
modes with long relaxation times r T are responsible 
for most of the dissipation and thus control the drift ve- 
locity. Steady-state motion has an infinite characteristic 
time T, so that only the zero mode - a rigid translation of 
the wall - with t = 00 is relevant in this case. As shown 
in Sec. IIII Bl the oscillatory regime has a characteristic 
time scale T — 7r/(7/io-ff). This regime has one addi- 
tional soft mode with ti > T as long as the field is not 
too strong. As the field strengthens, T becomes shorter 
and eventually additional modes become soft, T2 ^ T etc. 

A large gyrotropic force creates a softening effect. Con- 
sider, as an illustration, system ([5]) with two modes ^1^2, 
free energy U = +^2)/2, viscosity matrix Tij = TSij 
and gyrotropic matrix Gij = Gsij, where Sij is the an- 
tisymmetric tensor with £12 = -fl. When G = 0, one 
has two purely relaxational modes with r = T/k. These 
modes will be soft for small stiffness k or high viscos- 
ity r. In the opposite limit G ^ F the solution ex- 
hibits underdamped oscillations with the relaxation time 
T = GV(fcr). The latter exceeds the G = result by 
a factor (G/F)^ :§> 1. We will see that in permalloy, 
the smallness of a « 0.01 means that G ^ F, i.e., the 
gyrotropic force indeed dominates the viscous forces. 

The general formalism described in this section is illus- 
trated in Appendix on the classic problem of a Bloch 
domain wall first considered by Walker— and recently 
reviewed by Thiaville and Nakatanijii 



TABLE L Dimensionless ratios Tij /G for permalloy strips of 
widths w = 200 and 600 nm computed for the vortex wall 
model of Youk et al}^ with vortex core placed in the cen- 
ter of the strip. Changing the position of the vortex does 
not change the results drastically. The ratios depend on the 
vortex chirality x = ±1- 





Txx /G 


Txy/G 


Fyy/G 


w = 200 nm 


0.044 


0.031X 


0.049 


w = 600 nm 


0.116 


0.103X 


0.131 



III. TWO-COORDINATE APPROXIMATION: 
GENERIC FEATURES 

For domain walls under consideration, a large gy- 
rotropic force is associated with the motion of the vortex 
core (Appendix |B]). As a result, one may expect the 
vortex core motion to represent the softest modes of the 
system, which is indeed the case. We therefore consider a 
minimal description of the vortex wall with just two coor- 
dinates {X, Y) giving the location of the vortex core.^^ In 
order to calculate, e.g., the viscosity tensor Fy , we must 
have a model m — m{X{t),Y{t)) for the wall. Such a 
model is discussed in Appendix [Cl and the resulting vis- 
cosity components for permalloy strips of widths w = 200 
and 600 nm are tabulated in Table |T1 

However, once we have settled on the X and Y posi- 
tions of the vortex core as our collective coordinates, we 
may draw some general conclusions regarding the mo- 
tion that are independent of the model we choose for the 
wall. In particular, this choice of coordinates leads di- 
rectly to a universal time for the vortex to cross the strip 
in the transverse direction, in agreement with experimen- 
tal observations."^ 



A. The gyrotropic tensor 

Because we are using only the two coordinates of the 
vortex core to describe the wall motion, the antisym- 
metric tensor Gij has a single independent component 
GxY = —Gyx- As shown in Appendix [Bl as long as the 
vortex core is rigid, Gxy has a universal value pG, where 
G = 27r Jt, J = Mf-f is the density of angular momen- 
tum, t is the thickness of the film, and the polarization 
p = ±1 indicates the sign of the out-of-plane component 
of the magnetization within the vortex core. 



B. X-dependence of the free energy 

A vortex wall with the vortex core at {X, Y) has the 
free energy U ~ U{X,Y). Because of the translational 
invariance along the length of the strip, the dependence of 
the energy on the longitudinal coordinate A" at a fixed Y 
is trivial. A rigid shift of the wall by AX alters the length 
of the two domains with the opposite magnetizations by 
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+AX and —AX and thus changes the Zeeman energy 
-Aio/H ■ MdV by -QHAX, where Q = 2fj,oMtw is 
the magnetic charge of the domain wall (see Appendix 
O Fig. [TO] in). Hence, 



U{X,Y) = -QHX + U{Y). 



(9) 



Note that the longitudinal force Fx — —dU/dX — QH is 
independent of the transverse coordinate Y and in fact of 
any other details of the wall structure. This has interest- 
ing consequences for wall motion in very high magnetic 
fields. In this limit both gyrotropic and Zeeman forces 
are much larger than viscous ones, i.e., such a regime 
corresponds to the limit of zero viscosity. If viscous 
forces are completely neglected, the longitudinal Zeeman 
force must balance the longitudinal component of the gy- 
rotropic force pGY giving a constant transverse velocity 
of the core Y = —pQH/G — —p^^oHw/n and a univer- 
sal timei^iii for the vortex to cross the strip. 



T : 



w 

W\ 



(10) 



The transverse coordinate Y thus oscillates at the Lar- 
mor precession frequency lul = lHoH . It is a remarkable 
fact that in zero viscosity limit the same frequency is ob- 
tained for a completely different domain wall in Walker's 
problem (see Appendix 1X1 and Ref. ITsh. 



C. y-dependence of the free energy 

In the absence of the applied field the system possesses 
a symmetry with respect to 180°-rotations of the strip, 
X ^ -X SinAY ^ -Y, so that U{Y) w kY^/2 to 
the lowest order in Y. (Refiectional symmetry Y i-^ —Y 
is broken by the chiral nature of the vortex wall.) An 
applied magnetic field breaks the rotational symmetry 
allowing a term linear in Y . We thus have 

U{X, Y) = -QH{X + rxY) + kY^/2 + 0{Y^), (11) 

where x = ±1 is a chirality of the vortex and r is a numer- 
ical constant, and k characterizes the restoring force that 
tries to keep the vortex in the center of the strip. This 
results in the transverse force Fy ~ rxQH — kY. The 
model we use for the vortex wall shows a significant trans- 
verse Zeeman force with r of order 1 (see Appendix \G\ 
Fig. [TO)) , and nonzero restoring force. 



D. One-dimensional effective model 

In the absence of dissipation, the equations of motion 
may be derived from the Lagrangian 



L{X, Y) = -pGXY - U{X, Y). 



(12) 



Note that dL/dX = —pGY, so Y is the canonical mo- 
mentum conjugate to X. 



If we are interested only in the longitudinal motion 
of the vortex, we may eliminate the transverse coor- 
dinate Y using the equation of motion. In the har- 
monic approximation for the energy of the wall (jlip . 
Y = {rxQH - pGX)/k. Substituting this into Eq. ([T2|) 
gives the effective Lagrangian representing a massive par- 
ticle in one dimension subject to a constant external 
force, 



L{X) = 7nX'^/2 + QHX, 



(13) 



with the effective mass m — G^/k. Using our estimate 
of the stiffness k, Eq. (|C17p . and the definition of G 
in Sec. IIII Al one finds that this mass is typically of 
order 10^^^ kg. This is the same order of magnitude 
as that found experimentally for a transverse wall by 
Saitoh et alJ^ and theoretically for a one dimensional 
wallji2iS£ The mass increases with the width of the strip 
and changes only weakly with the thickness. 

The effective one-dimensional description shall prove 
useful in Sec. |Vl where additional degrees of freedom af- 
fect the motion of the wall. If these modes also have 
restoring forces acting on their conjugate momenta, the 
model becomes one of interacting massive particles in one 
dimension. 



IV. TWO-COORDINATE APPROXIMATION: 
DYNAMICS 

The equations of motion ([5|) for two generalized coor- 
dinates ^0 = X and £,i — Y read® 



(14) 



where Fy ~ Tji is the viscosity tensor, p is the polar- 
ization of the vortex core, and Ey is the 2x2 antisym- 
metric tensor with Eqi = +1. The generalized forces 
Fi — —dU/d^i are derived from the free energy pT|) . We 
thus arrive at the following equations of motion for the 
vortex core: 

+ ,15.) 



Y = 



Txx G2 + dctr 



G2 + det r 

where the equilibrium Y position of the vortex is 
{X9-P)GQH 



Y — 

^ eq 



fcF 



XX 



with 



{rVxx - X^xy)/G « 1 



(15b) 



(16) 



(17) 



It is worth noting that the magnitudes of the transverse 
displacement jVcql are slightly different for the two pos- 
sible values of the product px of the polarization and the 
chirality. This effect can be traced to the lack of the re- 
flection symmetry y i— > —y in a vortex wall, which leads 
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He and is guaranteed by the smallness of parameter 



Small dissipation regime 



Higti-field constant 
mobility regime 

c> ci> ci> 




FIG. 2: A sketch of the domain wall drift velocity V{H). 
Lower panels show the vortex motion trajectory (solid lines) 
and equipotential lines (dotted) at different field magnitudes. 
(A) - below the Walker breakdown, (B) - just above the 
Walker breakdown, (C) - small dissipation regime, (D) small 
dissipation regime with dominating Zeeman force. 



to nonzero transverse components of the Zeeman force 
rxQH and the viscous force TyxX. As a result, the tra- 
jectories of vortex cores with the same polarization and 
opposite chiralities x ^-re slightly different. 

Analysis of equations of motion (fT5|) yields three dis- 
tinct regimes (Fig. [J). Below a critical field He we find 
steady viscous motion with a high mobility fi — dV/dH 
(see Sec. IIV Al below) . Immediately above the critical 
field He the vortex motion becomes oscillatory, and the 
drift velocity of the wall quickly decreases as the applied 
field grows (see Sec. IIV El for details). At much higher 
fields, H ^ Hq = kw/2Q, the drift velocity rises linearly 
again but with a much lower mobility than at low fields 
(Sections IIV Cl aud IIV DI below). The separation of scales 



A. Low fields {H < H^) 

In a low applied field the wall exhibits simple viscous 
motion. The transverse coordinate of the vortex will 
asymptotically approach its equilibrium position Feq, as 
long as the latter is within the strip. The wall then moves 
rigidly with a steady longitudinal velocity V = QH/Txx 
giving the low-field mobility 



//LF = dV/dH = Q/T 



XX- 



(18) 



Using the calculated value of Txx (|C11|) for a strip of 
width w = 600 nm and Gilbert damping a ~ 0.01 yields 
/iLF = 29 m Oe~^, which is not too far from the value 
of 25 ms~^ Oe~^ measured by Beach et al^ 



B. Critical field {H = He) 



As demonstrated in Table [J and Eq. (jClip . the viscous 
force is small in comparison with the gyrotropic one in 
permalloy strips with a width below 1 /im. As a result, 
the equilibrium of a vortex in the transverse direction is 
set mostly by the balance of the transverse components 
of the gyrotropic force GV and the restoring force — fcl^q. 
The low-field regime ends when the equilibrium position 
of the vortex core is pushed beyond the strip edge, |Foq| > 
w/2, making the steady state unreachable. The critical 
point is reached when |l^q| = w/2, so that the critical 
velocity is 



Vc 



yielding the critical field 



kw 
2G' 



Hr 



^xxVc kwTxx 



Q 



2QG 



(19) 



(20) 



Beach et alA found a critical velocity of 80 m/s in a 
permalloy wire 600 nm wide. With the aid of Eq. (jC17[) 
for the stiffness constant fc, we find Vc = 188 m/s in 
such a wire. This is about twice as high as that ob- 
served experimentally. It is, however, much closer than 
that in Walker's one dimensional model of the wall, 1770 
m/s.^'^'^^ It is notable that our estimate lies between the 
experimental measurements and the critical velocity of 
256 m/s, found in micromagnetic simulations by Yang et 

The critical velocity grows logarithmically with the 
width of the strip, and nearly linearly with its thickness. 
It is easy to see that the latter result is valid beyond the 
model of a vortex wall adopted in this calculation. The 
two forces balancing each other at Y^q = w/2 scale dif- 
ferently with t. While the gyrotropic force is linear in 
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the restoring force comes from the magnetostatic en- 
ergy, which represents Coulomb-hke interaction of mag- 
netic charges with density 0{t), hence (the dipolar part 
of) the restoring force is quadratic in t up to logarithmic 
factors. This implies that Vc is roughly proportional to 
t. 



C. General remarks for high fields {H > He) 

Our numerical simulations in strips of thickness t ~ 20 
nm and width w = 200 nm indicate that, after the origi- 
nal vortex with a core polarization p is expelled from the 
strip, a new vortex is injected at the same location with 
the opposite polarization —p. The vortex thus moves be- 
tween the edges switching its core polarization each time 
it reaches an edge. 

Once the transverse coordinate of the vortex Y be- 
comes a dynamical variable, the motion acquires an 
entirely different character. As pointed out above, in 
permalloy strips the gyrotropic force GV dwarfs the vis- 
cous force rV. To zeroth order in Txx/G, the dynamics 
is purely conservative. Using the general form of the La- 
grangian p2p one can infer that the vortex core moves 
along equipotential lines U{X,Y) = const. Absent vis- 
cosity, the vortex would oscillate back and forth with the 
crossing time given by Eq. (fTO|) but the wall would not 
move on average. Any total x-displacement of the wall 
releases Zeeman energy, and thus requires viscous fric- 
tion to dissipate it. At a small finite viscosity the vortex 
trajectory slightly deviates from the equipotential lines, 
and the wall slowly drifts in the longitudinal direction 
(Fig.ElCD). 



D. Very high fields (H > Ho) 

We first demonstrate that at very high fields the ve- 
locity is again proportional to the field and calculate the 
high-field mobility for two-coordinate models. The new 
field scale Hq is set by the requirement that the restoring 
force —kV be negligible in comparison with the Zeeman 
force QH. The characteristic field is 

Ho^'^^^»H.. (21) 

When H ^ Ho, the dynamics is dominated by the Zee- 
man and gyrotropic forces, so that the vortex moves al- 
most along an equipotential line X + rxY = const. 

As a result of the drift with a velocity Vd, the Zeeman 
energy goes down on average at the rate QHVd- It is 
dissipated through heat generated at the rate 



W^TY ^Y^ ( -rx 1 )f 




The transverse velocity of the vortex core reflects the 
balance between the longitudinal components of the gy- 
rotropic and Zeeman forces: Y k, —pQH/G. We thus 



flnd the drift velocity Va = fJ-npH with the high-field 
mobility 

MHF = Q (ryy - 2rxTxY + r^^xx) /G^. (22) 

The high-field (HF) mobility ((22l) is suppressed in com- 
parison to the low-field (LF) one ([T8|) by a factor of 
0(r2/G2) = 0{a^). In the experiment of Beach et al.± 
fJ-HF / fJ-LF ~ 0.1. Using the F values for the more realistic 
model of Youk et al. (shown in Table |T]) and r « 2 as pre- 
dicted in the AppendixO we find /xhf / Mlf ~ 0.01. Since 
the predicted low-field mobility matches the experimental 
value fairly well (see Sec. lIV A)) , the calculated high-field 
mobility is much smaller than the observed value. 



E. High fields {H > He): Details 

To find the drift velocity of the vortex at fields above 
the vortex expulsion field we determine the total x- 
displacement of the vortex over a full cycle of motion from 
the top of the strip to the bottom and back again. The 
crossing time will be slightly different on the upward and 
downward trips due to a nonzero transverse component 
of the Zeeman force (Sec. IIII Cp . 

Solving Eq. with polarization p — ±1 gives us the 
crossing times and displacements Ar+ and AX^ (top to 
bottom) and Ar_ and AX^ (bottom to top): 



AX± 
AT± 



QHAT± G T TxY 



= 2- 



Fxx 
.G^ 



^xx 



-w, 



kr 



det r , Hc± 

arctanh— — — . 

H 



XX 



(23a) 
(23b) 



where we define Hc± = He/ {I T X5)j with g defined by 
Eq. (fT7|l . One can see that Eq. (|23bp reduces to the 



universal crossing time of Eq. 
The drift velocity is 



T0| in the hmit F, 



0. 



Vd = 



AX+ + AX- 
AT+ + Ar_ 



H 



2(l + detF/G2 



arctanh(iJc+/-H^) + a.Tcta.nh{Hc- / H) 



(24) 



Note that the critical fields Hc± are slightly different re- 
flecting a coupling between the vortex polarity p and chi- 
rality x seen in Eq. (|16p . 

An expansion of Eq. p4)) in powers of 1/H yields the 
high-field result ([22|l . The same expansion allows to find 



the field -ffmin Sit which the domain wall velocity has a 



mmimum: 



He 



V3(.g2 + detF/G2) 



(25) 



Using expressions (j23al) for the displacements AX± it 
is possible to characterize the fields at which the vortex 
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FIG. 3: (Color online) The drift velocity Vd of the domain 
wall as a function of the applied field H for a permalloy strip 
of width w — 200 nm and thickness t = 20 nm. Dashed 
vertical lines mark the critical fields He- and Hc+- 



trajectory approaches the zero-dissipation limit. Typical 
trajectories of the vortex are shown in Fig.[5]JC,D). They 
are close to the equipotential lines U{X, Y) = const when 
the displacement AX++AX^ in one back and forth cycle 
is negligible compared to the strip width w. Expansion 
in gives a criteria 



2G 



(g^ + detr/G^ 



2G 
STxx 



H 



< 1 



or 



2G 
H 



detr/G^) < 1 



> 



2G 



3r 



XX 



(26a) 
(26b) 



For r < 1 the first inequality is always satisfied when 
r < G holds. Eqs. ^ and ((26b)) show that vortex 
motion becomes approximately dissipationless for fields 
above the velocity minimum. 

In Fig. [3] the predicted drift velocity is compared with 
the results of numerical simulations for a permalloy strip 
of width w — 200 nm and thickness t = 20 nm.— 

The components of the viscosity tensor used are the 
predicted values for the Youk et al. model listed in 
Table HI The stiffness constant k of the restoring po- 
tential could not be calculated accurately because two 
of its main contributions, a positive magnetostatic term 
and a negative term due to Neel-wall tension, are of the 
same order of magnitude. This is not surprising given 
the proximity of the strip used in the simulations to a 
region where the vortex wall is unstable.— Instead, the 
relaxation time ri = (G^ -I- detr)/krxx was extracted 
directly from the numerics by fitting Y{t) to the expo- 
nential time dependence described by Eq. (|15p . i.e. 



Y{t) - Feq + [Y{0) ~ Feq] exp(-t/ri). 



(27) 



see Fig. For applied fields from 4 to 60 Oe, ti was in 
the range from 8.5 to 9 ns. This leads to a fc value roughly 
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FIG. 4: (Color online) The time constant n is determined 
by fitting the numerically calculated time dependence of the 
F-coordinate to a decaying exponential, see Eq. (|27p . 
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FIG. 5: (Color online) Numerical simulations above the crit- 
ical field show repeated oscillations in the width of the wall, 
measured as the distance between the half-antivortices. The 
patterns for various field values are offset vertically by 150 
nm for clarity. Each wall width is oscillating near w = 200 
nm. The vortex crossing time T and the decay time of wall 
width oscillations T2 are shown by arrows. 



twice that predicted by Eq. (IC17p . In accordance with 
Eq. (jlSp . the equilibrium position Y^^^ scales linearly with 
H. In calculating the critical velocity 14 — kw/{2G), it 
was necessary to replace w with an effective strip width 
WcS — w — 2R{) to account for a finite size of a vor- 
tex coreSi and edge defect. As the vortex approaches 
the edge, the short-range attractive exchange interac- 
tion overwhelms the dipolar potential and draws the vor- 
tex into the edge defect. This leads to a deviation of 
the transverse vortex coordinate plotted in Fig. 3] from 
Eq. (P7)) as the vortex comes close to the edge. From the 
observed trajectories we estimated i?o ~ 10 nm. 



V. BEYOND THE TWO-COORDINATE 
APPROXIMATION 

At low and intermediate fields, the two-coordinate 
model shows good agreement with numerical simulations. 
However, in higher fields the calculated drift velocity falls 
below the value observed in the numerics (Fig. [3]). The 
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discrepancy is caused by the softening of additional de- 
grees of freedom that have so far been ignored. The new 
soft modes provide additional channels for the dissipa- 
tion of Zeeman energy leading to an increase in the ve- 
locity. In fields above H2 — 35 Oe, the oscillation period 
T = n/{'-ffioH) becomes comparable to the decay times 
r of some of the harder modes. The new modes can be 
seen in oscillations of the total width of the wall that 
occur after the vortex has been absorbed by the edge 
and reemitted (Fig. [5]). Including additional modes in 
the calculation clarifies some of the significant features 
of the wall dynamics observed in numerical simulations 
at high fields. 

We first present a hypothetical toy model with two 
additional degrees of freedom in Sec. IV Al It provides 
a pedagogical example of coupled hard and soft modes 
with gyrotropic forces and sets the stage for a more real- 
istic model with six degrees of freedom appropriate for a 
vortex domain wall in higher fields, which is the subject 
of Sec. rVBl 



A. Pedagogical model with 4 coordinates 

The high-frequency oscillations of a vortex wall (Fig.O 
can be reproduced in a simple phenomenological model 
where the wall is treated as a vortex-antivortex pair with 
cores located at (xi, yi) and (x2, 112) ■ (A real vortex wall 
contains a vortex and two half-antivortices.^) The gy- 
rotropic tensor of the system has the following nonzero 
components: 

GxiVi = -Gy^xi = Gl, Gx2y2 = -Gy2X2 = ^2 . (28) 

Both defects are confined in the transverse direction 
and coupled to each other in the longitudinal direction 
with an equilibrium separation a. They carry magnetic 
charges Qi and Q2- The potential energy of the system 
is 



k{xi — X2 — a)^ 



The Lagrangian is 



-Gixiui - G2X2y2 - U. 



(29) 



(30) 



In this pedagogical example we assume that the viscous 
force acts on the antivortex only, so that the dissipa- 
tion tensor has a single nonzero component Tx2X2 = 
r > 0, giving rise to the Rayleigh dissipation function 
R = TX2/2. Also, it is important to note that the vortex 
core is flipped whenever it reaches the edge of the strip, 
yi = ±w/2, altering the sign of Gi. 



kiyi = —GiXi. The resulting dynamics is that of two 
massive particles moving in one dimension. The potential 
energy of the original model ([29|) translates into a sum of 
kinetic energy K and potential energy U of the massive 
particles: 



E = K + U 



mix\ ^ m2x\ ^ k{xi 



X2 



H{QiXi + Q2X2), 

(31) 

where the masses are determined by the gyrotropic and 
stiffness constants, nii = Gf/ki. Continuous evolution 
of the system is disrupted when the velocity of the vor- 
tex attains a critical magnitude, — Vc ~ kiw/2Gi, 
signaling that the vortex has reached the edge. At that 
moment Gi changes sign and the longitudinal vortex ve- 
locity ii = — fcij/i/Gi is reversed, xi — > — ii. After that 
continuous evolution resumes. Nothing special happens 
to the antivortex at that moment. 

The natural modes of the system are the center of mass 
X = (mixi + 11-12X2) /M and relative position a; — xi — 
X2 — a. Their energies are 



Ec 



MX' 



M 



2 
;.2 



QHX, 



lix" kx 
-Bid = — + — qHx, 



where M — mi + m2 and fi 



(32a) 

(32b) 
are the 

total and reduced masses and Q = Qi + Q2 is the to- 
tal magnetic charge. The value of the relative magnetic 
charge q will not matter to our analysis. 

We consider the limit in which the antivortex is 
strongly confined in the transverse direction, so that 
k2 ^ ki and as a result, the antivortex is much lighter 
than the vortex, 



(33) 



Here we rely on the fact that the gyrotropic coefficients of 
the vortex and the antivortex have the same magnitude, 
see Appendix [Bl 

The center-of-mass position X is a zero mode with an 
infinite relaxation time. Its velocity X is an overdamped 
mode with the relaxation time ri ^ toi/F = G\/{kiT), 
which is essentially the same as in our two-coordinate 
model of the vortex wall. The remaining two modes rep- 
resent underdamped oscillations of the relative coordi- 
nate and velocity with the frequency uj ~ and re- 
laxation time T2 ^ 2r7i2 /F <C ti . Thus x and x are much 
harder than X and X and there is a regime, ti ^ T 3> T2, 
where we may neglect the relative motion as a hard mode. 
In that regime, the dynamics of the wall reduces to that 
of our simple two-mode model described in Sec. IIVI 



Two massive particles 



2. Losses from the hard modes 



It is instructive to eliminate the transverse coordinates 
yi and y2, as in Sec. lIIIDl using the equations of motion, 



Let us compare the energy losses from the soft and hard 
modes to check whether our neglect of the hard modes is 
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justified. In the limit T2 <?C T, the relative position x and 
velocity x quickly reach their equilibrium values retaining 
them through most of the crossing period T. However, as 
the velocity of the vortex ±1 reverses its sign, the relative 
velocity x = xi — X2 changes as well increasing from zero 
to ±2wc- The energy of relative motion increases by 
at the expense of the center-of-mass energy. (Note that 
the total energy (pij) does not change at all since it is 
not sensitive to the sign of ii.) The kinetic energy of 
relative motion will be fully dissipated during the 

initial stage of the next crossing. 

The energy lost by the oscillatory relative motion, 
AE2 = 2^v^, should be compared to the total loss in- 
curred primarily through the overdamped main mode, 



AE = 



TX-'dt . 



TviT 



(34) 



The fraction of energy dissipated through the hard modes 
is thus 



AE2 
AE 



3t2 
T ■ 



(35) 



As long as relative motion remains hard, T2 <^ T, energy 
loss associated with it is negligible and the motion of the 
system is well approximated by the soft modes X and X 
alone. 



B. Model with 6 coordinates 

Transient oscillations of the hard modes found in the 
pedagogical model are similar to the oscillations of the 
width of the vortex wall (Fig. (5]), measured as the dis- 
tance between the two edge defects. In our numerical 
simulations, after the vortex is absorbed and reemit- 
ted at the edge, the width of the wall displays under- 
damped oscillations. In moderate fields, these oscilla- 
tions decay during the crossing time T and are reacti- 
vated at the beginning of the next cycle. This activation 
is evidence of energy transfer from the vortex core to 
harder modes, just as occurred in the pedagogical model 
of Sec. lV A[ In sufhciently high fields, the vortex crossing 
time T = 'k/{'^plqH) becomes comparable to the relax- 
ation time of the new modes. Notably, this occurs in the 
same region of field strength in which the two-coordinate 
model prediction for the drift velocity begins to deviate 
from the results of numerical simulations (Fig. [3]). New 
dissipation channels become important in higher fields, 
changing the dynamics of the wall. 

Further evidence of the need to include additional 
modes may be gleaned by comparing the equations of mo- 
tion of the two-coordinate model to the numerical simu- 
lations directly. The two-coordinate model predicts that 



pGY = QH-TxxX. 



(36) 



This equation expresses the change of linear momen- 
tum in the x direction due to external forces, and holds 
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FIG. 6: (Color online) Top row: The left hand side of Eq. ^ 
for the two-coordinate model plotted for the first three passes 
of the vortex across the strip under an applied field of 32 
Oe. Bottom Row: When a finite mass m = 1.1 x 10"^'^ kg 
is included for the edge defects, as in Eq. I|47|l . the total x- 
momentum of the modes is conserved. 



while the vortex is inside the strip. Indeed, the momen- 
tum canonically conjugate to coordinate X is dL/dX = 
—pGY p2)) . As a consequence of the translational invari- 
ance of the strip in the x direction, this momentum is in- 
fluenced only by forces external to the wall, i.e. the driv- 
ing force QH from the field and the drag force —TxxX. 
Equation ((36)) may be integrated to yield 

^X -pY-^t^ const (37) 

The above holds while the vortex is crossing the strip. 
Formally Eq. ([57]) does not describe the collision with 
the edge, but it can be checked that the effect of collision 
and corresponding vortex polarization flip is captured by 
a change of constant on the right hand side of Eq. ([57)1 
by const —^ const ± w. 

We can test this prediction using numerical simula- 
tions. The value of the left-hand-side of this equation is 
plotted in the top row of Fig. [B]for the first three passes of 
the vortex across a 200-nm-wide, 20-nm-thick strip when 
a 32-Oe driving field is applied. For Txx, we use the the- 
oretical estimate Txx = 0.044G for the model of Youk 
et al. (Table IJ). In the two-mode approximation, the 
plotted value should be constant during each pass. How- 
ever, it is clearly not constant and displays oscillations 
reminiscent of those in the wall width (Fig.[5|). 

Equation (|37p. derived in the two-mode approximation 
states that momentum along the nanostrip is influenced 
only by two external forces. Its violation reflects a trans- 
fer of momentum between the vortex and other modes 
of the wall. In a model with a larger number of modes 
a similar equation on the x-direction momentum can be 
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derived and checked numerically. It provides a good test 
showing whether the new set of modes is sufficient to 
capture the actual wall motion. 

The evidence in Fig.[5]and Fig.[n]shows that the modes 
associated with the positions of the edge defects become 
important to the motion of the wall in intermediate to 
high fields. We add the coordinates of these defects to our 
description in order to explain qualitatively their behav- 
ior as observed in numerical simulations. We use some 
simple modeling for the forces holding the edge defects in 
place relative to the vortex core. As with the vortex itself, 
each edge defect is described by two coordinates that are 
coupled to one another through the gyrotropic tensor. In 
the case of the vortex, these two variables are the X and 
Y coordinates of the core. For the edge defects, we use 
the X-positions X± of the defects and the out-of-plane 
angles 0± of the magnetization at their cores. It is not 
surprising that a gyrotropic coupling should be present 
between these coordinates, as it is the same pairing that 
occurs in Walker's problem of the Bloch wall — see Ap- 
pendix A for details. 

1. Energy 

In order to describe the motion of the two edge defects, 
two springs are added to the energy U{X, Y) in Eq. (fTTj) 
that act to keep the upper {X+) and lower edge 
defects on a 45° line with the vortex 

t/e. = |(x + xr-X,Tf)^ (38) 

Here again x is the chirality of the vortex. (Which of the 
two edge defects is ahead of the vortex depends on the 
vortex chirality.) 

To take into account the energy cost of bringing mag- 
netization out of the plane of the strip near the edge 
defects, the restoring springs are added for the angles 6*+ 
and 9- associated with the upper and lower edge defects, 
respectively: 

Ue± = ^el. (39) 



; ^^^^^^^^^^^^^^^^^^^^ 



are added to the Lagrangian. 

The gyrotropic coefficients for the upper and lower 
edge defects have opposite signs because they wind in 
opposite directions as one moves from left to right. This 
effect is visible in simulations even below the first break- 
down field (Fig. [7]) . As the wall moves to the right at a 
steady velocity, the lower defect tilts up out of the plane, 
while the upper defect tilts into the plane. 

The total Lagrangian for this extended model is then 

L = -pGXY - U{X, Y) + Y^(Lt,- Ueb - Ueb), (41) 
where U{X,Y) is given by Eq. (fTT]). 

2. Viscosity 

It remains to determine the components of the viscos- 
ity tensor for the three defects. Under the assumption 
that the three defects may move independently of one 
another, it is logical to assume that the viscosity tensor 
F will be essentially diagonal. Each of the edge defects 
carries with it a portion of the wall, and that portion's 
associated viscosity. Since the vortex is moving indepen- 
dently of the edge defects, there is no asymmetry that 
would lead to off-diagonal terms. The most general (sec- 
ond order) Rayleigh dissipation function consistent with 
these conditions is: 

R^lY.^.i^\ (42) 

i 

where Tx^ = = Te, and Tg^ = T0_ — Tg. The 

largest dissipation component is expected to come from 
the X± positions of the half-antivortices, as they influ- 
ence the motion of the Neel walls that emanate from them 
(see Fig. [TT|) . These Neel walls carry most of the viscosity 
of the wall, as they are the regions where the magneti- 
zation changes most rapidly with position. According to 
Eq. (j6bp , the viscous tensor components have the largest 
contribution from the regions where the change of magne- 
tization with the coordinate is the largest. The smallest 
dissipation component is expected to be that associated 
with the out-of-plane angles of the edge defects. These 
angles influence only a small region of magnetization im- 
mediately around the cores of the half-antivortices, and 
so will not have large contributions to the integral in 
Eq. (j6b[) . Therefore, we ignore Tg in what follows. 



FIG. 7: (Color online) Image of a wall moving to the right 
with an applied field of 4 Oe. Shaded regions indicate out-of- 
plane magnetization. 

To represent the gyrotropic coupling of the positions 
of the edge defects to their out-of-plane angles, the terms 

L± = TG2X±e± (40) 



3. Connection with the two-mode approximation 

This model can be related back to the model in which 
only the vortex is free to move by letting fc2 — *■ oo and 
kg oo. This wiU require X++X^ = X+xY, -x{X+- 
X_) = w, and 9^ = 9_ = for all time, leading to the 
reduced Lagrangian and Rayleigh function for just the 
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FIG. 8: (Color online) Top: Numerically determined x- 
positions of the vortex and the two half-antivortices in a 20 
nm thick, 200 nm wide permalloy strip as functions of time for 
H = 32 Oe. Bottom: The longitudinal positions of the vortex 
and the two half-antivortices in the model with massive edge 
defects. The vortex merges with an edge defect each time it 
hits the edge. Note that the oscillation of that edge defect is 
suppressed compared to the other. 



vortex position: 



(43a) 



i?, = VxX'l2 + VYY'l2 + V,[X + xy) ■(43b) 

Identifying Vx = ^xx - ^xy, Ty = Vyy - ^xy, 
and Fe = xF^f/Z, these are exactly the functions used 
when only the components of the vortex position are soft 
modes (Sec. IIVI) . It becomes clear then that the off- 
diagonal elements of the viscosity tensor in the 2-mode 
model arise directly from the asymmetric equilibrium po- 
sitions of the edge defects. By comparing the identifica- 
tions above with the viscosity values in Table [T] and Ap- 
pendix [Cl one can see that Fg is indeed the dominant 
viscosity term, especially as the width increases. Now 
that we have determined the viscosity values for the new 
model in terms of the old one, we can check the descrip- 
tive value of the model by testing the conservation of the 



total momentum in the numerical data. This momentum 
is derived in the next section. 



J^. Mass and momentum of the edge defects 

If we care only about the x-positions of the edge 
defects, we may integrate out the values of 9± as in 
Sec. nil Dl to produce masses m = G\/ke for the edge 
defects. The parameters kg and G2 then appear only in 
this combination. In effect, these masses occur because 
the moving wall is able to store kinetic energy in the out- 
of-plane angle of the edge defects. The Lagrangian then 
becomes: 

L = -pGXY - U{X,Y) + Y,{mXll2 - Ueh)- (44) 

The total x-direction momentum is then equal to 

= -pGY + m{X+ + X^) 

It is influenced by two forces in the x-direction. First, 
viscous drag acts on each of the defects, totaling to 



Fd — —TxX — TeX+ — TeX^ 



(45) 



Second, the external field applies a constant force Fh = 
QH on the wall. 

As in the two coordinate model, the translational in- 
variance of the system implies that the internal forces 
between the defects will have no influence over the total 
X momentum. Note in particular that the analysis here 
is independent of the spring constants k and k2- Thus 



dt 



= Fd + FH, 



This equation may be integrated to yield 

+ VxX + Fe (A+ + A:_) - QHt = const. 



(46) 



(47) 



We may determine the mass of the edge defects by fitting 
the left-hand side of Eq. P7)) to a constant. The left-hand 
side of this equation divided by the gyrotropic constant G 
is plotted in the bottom row of Fig. [6l The figure shows 
that Eq. gZl) with rn = 1.1 x lO^^^kg is satisfied by the 
numerical simulation at 32 Oe to nearly the accuracy of 
our simulation, which had a 2 nm cell size. 



5. Oscillations of the wall width 

Insight into the width oscillations can be gained imme- 
diately from the form of L and R in Eqs. ((4T|) and (|42p . 
Defining the wall width d by the difference — x(Ar+ — A_) 
between the forward and rear edge defects reveals that 
the conjugate pair (d, 9), where 9 = (0+ -I- 9-)/2, decou- 
ples from the remainder of the motion. Any deviation 
in d from w will decay over time with a time constant 
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FIG. 9: (Color online) Oscillations in the width of the wall in 
a strip with thickness t = 20 nm and width w = 200 nm in a 
field of H = 32 Oe. The upper curve shows the prediction of 
the model using the same parameters as in Fig. [S] It is offset 
vertically by 100 nm from the numerically determined result 
for clarity. 



Td = 2m /Te, while oscillating due to the gyrotropic cou- 
pling with the average out-of-plane angle of the defects. 

The equations of motion associated with the system 
described above can be solved to find the trajectories 
of the three defects. In order to take the collisions of 
the vortex with the edges of the strip into account, the 
polarization of the vortex is flipped every time it hits the 
edge. The effective width WeS ^ w — 2Rq, with Rg — 
10 nm, is used to account for the short-range exchange 
interaction of the vortex with the edge. The polarization 
flip changes the sign of the momentum of the vortex, as 
in Sec. IV AT] In addition, when the vortex hits an edge, 
it actually must collide with one of the edge defects to be 
absorbed. The vortex is reemitted from that defect with 
the opposite polarization. As a result, the polarization of 
the edge defect changes as well. This causes oscillations 
in the wall width d not because d itself has changed, but 
because its conjugate momentum 9 has been altered. 

The solution shown in Fig. [5] and Fig. [5] uses a crude 
method to demonstrate this effect, simply reversing the 
value of the out-of-plane angle of an edge defect whenever 
a vortex hits it. Thus the vortex and the edge defect that 
it hits both have their momentum change sign during 
the collision. This is similar to the effect discussed in 
Sec. IV A 11 in that it causes a sudden transfer of energy 
from the main mode, associated with the vortex motion, 
to ancillary modes associated with the oscillation of the 
edge defects. This effect is seen in Figs. [5]and[9l as the 
wall oscillation is reactivated every time the vortex hits 
the edge of the strip. 

The X coordinates of the three topological defects over 
time are plotted in Fig. [H] for the numerical solution and 
for the model in this section. The parameters used for the 
model are: G^/{kTxx) = 8.75 ns as measured in Fig. 21 



m = 1.1 X 10~^^ kg as determined in Fig. [6l r = 1.1 
to match the points where the vortex hits the wall, and 
k2/k — 1.3 to match the frequency of oscillations. 

The oscillations in the positions of the edge defects in 
the model and numerical solutions are comparable. In 
particular, the model correctly captures the fact that the 
oscillations in the farther edge defect from the vortex are 
larger in amplitude than those of the nearer one. This 
comes from the combination of flipping the sign of the 
vortex momentum with the resetting of the angle of the 
impacted edge defect. The oscillations caused by each 
effect interfere with one another. This interference is 
destructive for the near defect and constructive for the 
far one. 

Figure [9] shows the oscillations in the width of the wall. 
As expected, these oscillations have the largest amplitude 
immediately after the vortex is emitted from one of the 
walls. The time constant of the decay is of the same order 
in both the model and numerical solutions. This agree- 
ment for the decay constant Td = 2m/Te is expected, 
since this constant depends only on the edge defect mass 
measured in Fig. [6] and the viscosity coefficient Fe, which 
is known from the two-coordinate model. 

However, Fig. [9] also reveals the limitations of the 
model described in this section. While this model is use- 
ful in explaining the oscillations of edge defects, and the 
choice of coordinates shows excellent agreement with the 
x-momentum equation (|47p. as shown in Fig[6l there is 
still a significant discrepancy between the the theoret- 
ical prediction and the numerical simulation in Fig. [S] 
The average wall width seen in numerical simulations 
changes as the vortex crosses the strip, an effect that 
is not predicted by the theory thus far. To account 
for this effect, additional terms must be added to La- 
grangian (|4ip . In particular, this model thus far has ig- 
nored field-dependent contributions to the interactions 
between the defects, as well as terms of higher than 
quadratic order in the coordinates. 



VI. DISCUSSION 

We have explored the dynamics of a vortex domain 
wall in a magnetic strip of submicron width. We have 
applied the method of collective coordinates^ to the case 
when the wall has two soft modes related to the mo- 
tion of the vortex core. A simplified model of the vortex 
domain wall described in this paper yields solvable equa- 
tions of motion. The calculated mobility of the wall in the 
steady-state viscous regime at low fields agrees well with 
the value measured by Beach et alfi The steady motion 
breaks down when the equilibrium position of the vortex 
moves beyond the edge of the strip. The critical veloc- 
ity depends strongly on the magnetization length and the 
sample thickness, and weakly on the width; its calculated 
value agrees reasonably well with the data of Beach et alA 
The dynamics above the breakdown changes its charac- 
ter from overdamped to underdamped. In this regime the 
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velocity sharply declines at first but later starts to rise 
again as the field strength increases. The predicted high- 
field mobility is reduced in comparison with the low-field 
value by the factor //hf/a'lf ~ 0.01; the experimentally 
observed reduction is not as strong: /inp/MLF ~ 0.1.^ 

We have compared the results of our theory to numer- 
ical simulations of a moving vortex wall in a permalloy 
strip of width w — 200 nm and thickness t ~2Q nm. We 
have found that the results obtained via the collective- 
coordinate method are in quantitative agreement with 
the simulations both in the viscous regime at low applied 
fields and in the oscillatory regime at moderately high 
fields. As the field strength increases further, the veloc- 
ity observed in numerical simulations diverges from the 
predictions of the two-coordinate model. We have traced 
the breakdown of the two-mode theory to the appear- 
ance of four additional soft modes associated with the 
coordinates and out-of-plane angles of the edge defects 
(antihalfvortices) . 

As expected, the dynamics of a domain wall reduces 
to the motion of the elementary topological defects that 
make up the wall — bulk vortices and halfvortices at the 
edge. Of particular importance is the special kinemat- 
ics of vortices caused by their nonzero skyrmion num- 
bers (±1/2). On the formal level, the two coordinates of 
a vortex core form a canonically conjugated pair. As a 
consequence, when the vortex is constrained in the trans- 
verse direction by a parabolic potential, its longitudinal 
motion acquires inertia. By our estimates, in a strip of 
width w — 200 nm and thickness t = 20 nm a vortex 
wall has an effective mass of approximately 10^^^ kg. 
The mass of an edge defect is approximately 10"^'^ kg. 
Saitoh et al}^ determined the mass of a transverse wall 
in 70nm wide, 45nm thick permalloy to be 6.6 x lO^^'^ 
kg. 

In our analysis, we have made some simplifying as- 
sumptions that require further investigation. First, we 
have assumed that any vortex absorbed by the edge is 
immediately reemitted. While this holds true for our sim- 
ulated strips oi w = 200 and t = 20 nm, it is not always 
the case. There may be short delays between absorption 
and reemission during which the wall motion is again vis- 
cous; the higher mobility during this period would tend 
to increase the drift velocity. This type of motion has 
been observed in simulations of thinner strips by Lee et 

Second, while we have described the onset of wall width 
oscillations in Sec.|V]as a typical new mode, decaying os- 
cillations are not the only type of new mode that can 
occur. In particular, the number and dynamic charac- 
teristics of soft modes may change discontinuously as ad- 
ditional vortices or antivortices are created in the bulk 
of the stripf22, We have observed numerically the cre- 
ation and subsequent annihilation of a vortex-antivortex 
pair adjacent to the original vortex of the wall. As in 
the process described by Van Waeyenberge et al.r^ the 
pair creation mediates the flipping of the polarization of 
the vortex and results in the reversal of the gyrotropic 



force. Once the pair is created, the antivortex moves 
closer to the original vortex and together they form a 
topologically non-trivial bound state (a skyrmion) of the 
type described in detail by Komineasr& eventually de- 
caying in a spin wave explosion via a quasi-continuous 
process in which the conservation of topological charge 
is violated. ^^'^^ The new vortex, which has polarization 
opposite that of the original, takes over as the central 
vortex of the domain wall. This type of behavior may 
dominate the motion in very high fields, causing the pe- 
riodicity of the motion to deviate from the prediction of 

Eq. (nni). 
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APPENDIX A: WALKER'S SOLUTION FOR A 
180° BLOCK DOMAIN WALL 

We apply the method of collective coordinates to derive 
Walker's classic result^'^ for the dynamics of a 180° Bloch 
domain wall in a uniform magnetic field. The easy axis 
is z and the magnetization varies along the x direction, 
M = M.{x,t). The energy density per unit area in the 
j/z-plane is 

U = j dx[- ^JioHM cos e + KqSyu^ Bcos^ (j) 

-Kcos^e + A{e'^ + s\n^e(t)'^)], (Al) 

where primes denote derivatives with respect to a;, M = 
M(sin 9 cos 0, sin 9 sin (/), cos 9), H is an external magnetic 
field applied along the easy axis, K is the easy- axis 
anisotropy, and Kq — is the shape anisotropy 

refiecting the energy of a magnetic field — —M^ gen- 
erated inside the wall. 

On the basis of exact results for the steady state and 
numerical simulations beyond the steady state. Walker 
parameterized a domain wall by three collective coordi- 
nates: the center of mass A, the width of the wall A, and 
the (uniform in space) azimuthal angle (j). His Ansatz, 

cos0 = tanh(^^^S) , = 0(i), (A2) 

substituted into Landau-Lifshitz-Gilbert equations, 
yielded an approximate solution that reproduced the 
numerical results remarkably well. In Walker's solution 
the center of the wall X and the angle (p were treated 
as independent dynamical variables, while the width of 
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the wall A adjusted adiabatically to their instantaneous 
values. We will use the general formalism of Sec. |TI] 
to derive the equations of motion for the three modes 
and to deduce that the wall width is the hardest mode 
of the three with a parametrically short relaxation 
time, thus justifying the applicability of the two-mode 
approximation in the oscillatory regime where the period 
T and the relaxation times of the X, (f> and A modes 
satisfy the inequality 



oo. 



(A3) 



The equations of motion ([5]) require the calculation 
of generalized forces Fi and the matrices and Gij 



using the Ansatz for the shape of the domain wall (jA2[) . 
The antisymmetric gyrotropic matrix Gij has only two 
nonzero components, 



G. 



<j>X 



Gx<h — ±2 J 



(A4) 



for a domain wall with the asymptotic behavior cos ~ 
±1 at X — — oo. The viscosity matrix Tij is diagonal, 
with 



XX 



2aJ 



T^s = 2aJA, Faa = 



6A ■ 



Lastly, the generalized forces are 

Fx = 2/ioA/F, 

F^ = KQ/\s\n2(j), 

Fa = 2{A/A^ -K -Ko cos^ i 



(A5) 



(A6a) 
(A6b) 
(A6c) 



The width of the wall adjusts to the equilibrium 
value A(0) = AI[K + ifocos^ 0) on the time scale 
TA determined by its viscosity Faa and stiffness = 
-dF^/dA = 4A/A3: 



TA = 



Faa aiT^J 

^ 2A{K + cos2 , 



(A7) 



On longer time scales, the dynamical variables are the 
wall position X and the angle 0: 



Fx - TxxX + Gx4,'t> - 0, 



(A8a) 
(A8b) 



The wall position AT is a zero mode with an infinite re- 
laxation time. The terminal velocity is A = Fx /^xx = 
"/fioFlA/a. The angle reaches its equilibrium value, de- 
termined by the condition F^{(j)) -f G^xX = 0, on the 
time scale r^. Unlike the relaxation time of the wall 
width (jA7|l . is not determined by the angle viscosity 
F^^ and stiffness, = —dF^/dcj) — — 2AroA cos 2(/) > 0, 
alone. The gyrotropic coupling to A considerably softens 
this mode extending the relaxation time, 



The ratio of the relaxation times. 



T0 



Q;'^7r'^i4rocos2(?!) 



[l + oP-){K + KaCOii'^i 



< 



(AlO) 



is very small, 4 x 10 in yttrium iron garnet, where 
a « 10"^ and Kq/K = 2lJ^ It means that Walker's 
solution has a considerable range of fields H for which 
the period of oscillations T satisfies inequality (|A3ll and 
thus the system has two soft modes, A and (j). 



APPENDIX B: THE GYROTROPIC TENSOR 
FOR A BULK DEFECT 

The gyrotropic tensor is given by Eq. ([6]). If one de- 
scribes the motion of a bulk topological defect, such as a 
vortex or an antivortex with polarization p, it is conve- 
nient to choose the A and Y coordinates of that defect 
as the collective coordinates. This leads to a non-zero 
contribution to G, 



G 



XY 



J dVm- 



dm dm 

'dX^W 



(Bl) 



Note that the integrand here is closely related to the 
skyrmion density^ 



47r 



m • {dxm X dym) 



1 9(cos0, (/) 
47r d{x,y) 



(B2) 



In fact, if one makes the change of variables in Eq. (|Bip 
such that ^ ~* a-iid 'W ^ ~'5y^ integrand 

becomes A'Kp. 

This substitution is allowed only if the magnetization 
texture moves rigidly with the defect core. While this is 
not true everywhere within a domain wall, it is valid near 
the core itself, where the rigidity is forced by circulation 
of the magnetization. On these short distances, the ex- 
change interaction dominates and forces the magnetiza- 
tion to be out-of-plane, overcoming the shape anisotropy 
that would keep it in-plane. Because the exchange in- 
teraction is dominant, the shape of the core is mainly 
independent of its position. Fortunately, the rigid region 
near the core is the only region that contributes to the 
integral in Eq. (jBl|l . Because the integrand is a triple 
product, all three components of the magnetization need 
to be non-zero for a point to contribute. This only occurs 
near the core. 

Thus we find that Gxy = 47r J J pdV = 2nJtnp, where 
n is the in-plane (0(2)) winding number, which is 1 for a 
vortex and —1 for an antivortex, and p = Mz/\Mz\ = ±1 
is the polarization of the core. 



APPENDIX C: MODELING THE VORTEX WALL 



^2 

'-'<t,x 



Fj,rf,F 



XX 



2J(1 



kjT 



> 



2J 



XX 



uKq cos 20 C(Kq 



(A9) 



In this appendix we describe two models for the vor- 
tex domain wall based on two collective coordinates, the 
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^^'^'^'^'^'^'^'^'^'^'^'^'^'^'^'^'^'^'^'^'^'^'^ its simplicity, the model captures some of the essential 

^ s / features of the vortex wall, including the chiral properties 

^ ^ ^ ^ i / — that lead to off-diagonal terms in the viscous tensor and 

I / the transverse component of the Zeeman force on the 

vortex. 
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FIG. 10; A very simple (4-Neel-waIl) model of the vortex 
domain wall~ The vortex core is denoted by the filled circle. 
Four Neel walls separate regions of uniform magnetization. 
Shaded areas indicate a buildup of magnetic charge. The 
panels show states with a fixed longitudinal coordinate of the 
vortex X = const; the transverse coordinate is y = w/2, w/4, 
0, —wj'^, and —w/2 respectively. 



position of the vortex core (X, F). The simpler model 
provides a pedagogical introduction and the more sophis- 
ticated version is used to derive numerical values for the 
parameters used in Sec. IIVI 

In strips that support vortex domain walls, the domi- 
nant contribution to the energy in the absence of an ap- 
plied field is due to magnetostatic interactions. For any 
domain wall in a strip of width w, thickness i, and satu- 
ration magnetization M, there is a total magnetic charge 
Q — 2iJ,oMtw associated with the wall.^^ In a vortex 
wall, nearly all of this charge is expelled to the edge of 
the strip. We thus focus on models in which there is no 
bound magnetic charge in the bulk of the strip. 

We begin with a simple model in which a vortex wall 
consists of four straight Neel walls separating regions of 
uniform magnetization (Fig. llOp . Two of these walls cross 
at the vortex core, and each edge defect has two walls 
emanating from it at 45° angles from the edge. Despite 



FIG. 11: 3-Neel-wall model of the vortex domain wall pro- 
posed in Ref. I 111 . Solid lines denote straight Neel wall regions 
and dashed lines denote regions in which the Neel walls are 
parabolic. 



While the 4-Neel-wall model described above correctly 
determines the form of the free energy and the sign of the 
off-diagonal viscosity term Txy, it significantly overesti- 
mates the magnitude of the viscosity coefficients. This is 
mainly due to the inclusion of two Neel walls intersect- 
ing at the vortex core. In the more realistic model first 
proposed by Youk et al}^ one of these walls is replaced 
with two regions where magnetization rotates gradually 
(Fig. llip . In these curling regions, the magnetization an- 
gle is given hy (fi = /3-|-x7r/2, where P is the angular coor- 
dinate around the vortex core, originating at the x-axis, 
and X = ±1 is the chirality of the vortex. As a result, 
we are left with three Neel walls, with only one passing 
through the vortex. These walls are straight where they 
separate two uniform regions and parabolic where they 
separate a uniform and a curling region. This prevents 
bound magnetic charge in the bulk of the system. 

We proceed in this appendix to derive the components 
of the viscosity tensor and the values of the param- 
eters r and k used in the energy pT|) . In each case we 
begin with a brief explanation using the simpler model 
before deriving the numerical values using the model of 
Youk et al. 
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1. Viscosity tensor 

As long as the magnetization vector lies in the plane 
of the strip the definition of the viscosity tensor © may 
be expressed in terms of the azimuthal angle charac- 
terizing the magnetization: 



aJt 



2 dcj) dcj) 



(CI) 



where J — M/7 is the density of angular momentmn and 
t is the thickness of the film. 

The largest contribution to the viscosity comes from 
regions where the magnetization angle depends strongly 
on the collective coordinates X and Y. Note that 
the translational symmetry of the problem means that 
(j){x, X,y,Y) — (p{x — X,y,Y), allowing us to make the 
replacement d4>/dX — —d(f>/dx in all cases, regardless of 
the model. In the simple model (Fig. [TO|) an infinitesimal 
shift in X or y affects magnetization only in the vicinity 
of the Neel walls. 

We begin by considering the contribution to the vis- 
cous tensor of the Neel walls that emanate from the 
vortex core at ±45°. Near these walls, (f>{x,y, X,Y) = 
f{x — X ^ y ± Y). In these regions, derivatives with 
respect to collective coordinates may be reduced to ordi- 
nary gradients: d(j>/dX — —d(t>/dx — —f and d(j>/dY — 
—d<j>/dy — ±/'. As a result, the tensor components are 
equal to each other, up to a sign: 



i XX = 1^ FY Tl = aJt 



d xf 



(C2) 



Opposite signs of the off-diagonal component Txy can be 
understood by noting that, as the vortex moves along Y, 
the two Neel walls that pass through it shift along X in 
opposite directions, creating equal and opposite viscous 
forces in the X direction. 

As shown in Fig. 1101 the two peripheral Neel walls move 
in the same direction as the central Neel wall perpendic- 
ular to them. Since they are also of the same length as 
that wall, they have the same contribution to the viscous 
tensor. Adding the contributions of all four Neel walls 
yields a total 



XX 



r 



YY 



2xrxY, 



(C3) 



independent of the vortex position (for the simple 
model) . 

To obtain the absolute values, note that the expres- 
sions for viscosities (|C2p coincide, up to a constant fac- 
tor, with the exchange energy of the Neel wall, which has 
been calculated, e.g., in Ref. [sO- 

This energy may be determined by integrating the sur- 
face tension of the Neel wall, given by^S 



a^2A 



ov 



2V2A 



(sin (po - 4>o cos(/io) 



(C4) 



where z; is a coordinate perpendicular to the wall, and 
2(/)o is the angle of rotation across the wall. For a straight 
wall this angle is constant along the entire wall and the 
remaining integration along the wall is trivial. After an 
appropriate normalization we obtain viscosity coefficients 



^xx = Tyy = 2xTxY 0.608aJtw/A 



(C5) 



for the simple model, where the exchange length A — 
\/ A/ iiqM'^ =3.8 nm in permalloy. 

For a more accurate computation of the viscosity ten- 
sor, we turn to the model of Youk et ali^ (see Fig. [TT|). 
In this model, there are seven areas that contribute to 
the viscosity: 3 straight Neel wall segments, 2 parabolic 
Neel wall segments, and 2 regions of the bulk in which 
the magnetization curls. 

Eq. (jC2|) still holds for the straight Neel wall regions. 
However, one of the Neel walls passing through the vortex 
is now absent, and the outer walls have half the length of 
the inner. Their viscosity contribution is equal to that of 
the inner wall. As a result, the total contribution from 
the straight segments of Neel wall is 



r 



XX 



YY 



xr'xY 



0.30AaJtw/X. 



(C6) 



In a similar fashion, we may calculate the viscosity 
of the parabolic Neel wall segments. There is a slight 
complication in this case as these segments deform as 
the vortex moves across the strip. 

The parabolic Neel walls in Fig. [11] are described by 
the equation {x — X)'^ = (2y ± w){2Y ± w) for the up- 
per and lower regions, respectively.^^ Note that when the 
vortex moves straight up, the distance that a given point 
on the parabolic Neel wall moves is dependent on its y 
coordinate. For an infinitesimal displacement dY of the 
vortex upward from the center, a point on the parabolic 
Neel wall is moved to the left by dYy^l ± 2y/w. Thus, 
along these walls dcji/dY = -{d(l)/dX)^l±2y/w. 

In the case of the parabolic Neel walls, the angle 0o 
changes along the length of th e wall. For the lower wall, 
00 = (arctan(y/-y/w(u' — 2y)))/2 + 7r/4. This is equal 
to the angle of the wall normal away from the x-axis. 
Therefore, 



dx 



90 
dx 



V2 
A 



(sin 00 — 00 cos 00 ) cos 0o (C7) 



in the region near the lower parabolic Neel wall. We 
assume that 0o varies slowly on the scale of the Neel 
wall width. The variation of 0o along the wall leads to 
parametrically small, width independent corrections to 
the Neel wall viscosity calculated here. 

Making the usual replacement d(t>/dX = —d(j)/dx, we 
can perform the integration in Eq. (|C1[) numerically for 
each tensor component. When the vortex is centered on 
the strip, the contribution of the parabolic walls to the 
viscosity is given by 



= QmiaJtw/\, 
T^Y = 0.083aJiu;/A, 
T^Y = O.QmaxJtw/X. 



(C8a) 
(C8b) 
(C8c) 
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The final contribution to the viscosity comes from 
those bulk regions in which the magnetization curls. In 
these regions, it is most convenient to carry out the in- 
tegration of Eq. (jCl[) in cylindrical coordinates x — X ~ 
R cos P, y — Y — Rsinfi. When the vortex is centered on 
the strip, we find 



r„ 



aJt / dp / 

J-n/2 Jro 9^3 



-aJt / df] 



where tan /?o = — 1/\/8 and tq is the radius of the vortex 
core. 

The magnetization in these regions is given simply 
by = /? + 7rx/2, so that d(j)/dY = -(cos/3)/i?, and 
d(t>/dX = {sin 13) /R. Integrating (|C9[) numerically gives 



r^_y = aJt{^\n — - 0.398 



YY 



^ XY 



aJt ( - In 0.010 

4 To 



1 w 

= -axJt I - In 0.133 

2 To 



(ClOa) 
(ClOb) 
(ClOc) 



The symmetry of the domain wall when the vortex is 
centered on the strip implies that the contribution from 
both bulk curling regions is the same. Summing up all 
the contributions, we find 

Txx = aJd 0.418- + -In 0.797 , (Clla) 

V A 2 To / 

/ W IT W \ 

Tyy = aJt 0.470- + - In 0.020 , (Cllb) 

V A 2 ro / 

VxY = axJt (^0.442^ - In ^ + 0.265^ . (Cllc) 

It is instructive to compute the ratio of the viscous and 
gyrotropic forces Tij/G where G = 2-Kjt is the gyrotropic 
constant. Taking the vortex core radius ro = A = 3.8 
nm, we obtain the dimensionless ratios listed in Table [D 
The small value of Gilbert's damping in permalloy, a ~ 
O.Olfi^ leads to the dominance of the gyrotropic force in 
strips with submicron widths. The smallness of Tij/G 
can be exploited to organize an expansion in powers of 
this small parameter. 



2. Free energy 

In our simple model, magnetic charges form two lines 
of lengths w — 2Y and w + 2Y with constant charge den- 
sity per unit length p = fXoMt. As the vortex moves 
off-center, charge builds up on one edge of the strip or 
the other, leading to an increase in magnetostatic energy. 
The energy E{Y) = E{0) + kY^ /2 + 0{Y^) has a mini- 
mum at F = 0. (Here we chose F = to be in the middle 



of the strip.) This leads to a force —kY that acts to keep 
the vortex centered on the strip. 

As can be seen from Fig. 1101 transverse motion of the 
vortex core changes the total magnetization of the 
strip and thus affects its Zeeman energy. As the vortex 
core crosses from the top to the bottom edge of the strip, 
the Zeeman energy decreases linearly with the y-position 
of the core by a total of ApoHMtw^ . The dependence of 
the energy of the magnetic configuration on the vortex 
core position {X, Y) is therefore given by 



U{X, Y) = -Qi7(X + 2xY) + E{Y). 



(C12) 



The same principles hold true in the more realistic 3- 
Neel-wall model. As the vortex moves off center, the 
buildup of charge on one side of the strip causes a restor- 
ing force that pushes the vortex back toward the middle 
of the strip. As in 4-Neel-wall model, the Zeeman energy 
is given by —QH{X + rxY). Because the states of the 
system when the vortex is at either edge are the same as 
they are in the simplified model (see Figs. [TUl and [TT|l . 
the average value of r as the vortex crosses the strip is 
again 2. 

This is most likely an overestimate, as the states shown 
in Figs. [TO] and [11] when the vortex is at an edge have 
an extended edge defect with a core radius equal to the 
width of the strip. In reality, the half-vortex core is less 
extended due to the exchange cost of the surrounding 
Neel walls. This smaller core for the edge defect means 
that the total charge on the wall moves less in the X 
direction when the vortex crosses the strip than it does 
in cither of the models listed here. This leads to a smaller 
value for r. By fitting the vortex trajectory in numerical 
simulations to that of our collective coordinate analysis, 
we found r « 1.1 in a strip with the width w = 200 nm 
and thickness t = 2Q nm. 

In what follows, we derive the energy of the vortex do- 
main wall in the model proposed by Youk et alM- to the 
second order in the displacement of the vortex from the 
center of the strip. As a result, we determine the value 
of the stiffness constant k. There are four contributions 
to the energy in this model: the Zeeman energy of the 
magnetization in the external field, the magnetostatic en- 
ergy of the bound charges on the edges of the strip, the 
exchange energy of the bulk curling regions around the 
vortex, and the integrated Neel wall tension of the three 
Neel walls emanating from the topological defects. 

If (0, Y) is the position of the vortex, with Y mea- 
sured from the center of the strip, then the magnetiza- 
tion angles in the upper and lower curling regions are 
given by (/)± = arctan((y — Y)/x) =F 7r/2. Thus there 
are charge densities p±{x) = p^M sui(j)±{^w /2,x) = 
poM\x\/ x"^ + {w/2 =F F)2 spread out along the edges. 
The charge distributions are bounded in the center by 
the X position of the vortex and outside by the points 
L± = ^y/2w{w ± 2Y) at which the parabolic Neel wall 
segments hit the edges. The self-energies of these charged 
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FIG. 12; Value of the scaling function / for various values of 
the aspect ratio w/t. 



regions are given by 

SttiioJo Jo Jo Jo ^(x - .t')^ + (z - z')^ ' 

The interaction energy of the two charged segments is 
given by 



1 



dz I dz' I dx / dx' — -p 

47rpio Jo Jo Jo Jo V (a; + x'f + + {z - z'f 



Expanding to the second order in Y and integrating nu- 
merically, we find that the magnetostatic energy Em = 
Em{0) + ^V^. The value of /cm depends logarithmi- 
cally on the aspect ratio w/t. For w/t = 10, as in our 
numerical simulations, 



1.08 



w 



In general. 



Cm 



HoM-'t 



2i2 



(C13) 



where the value of / is plotted in Fig. [121 

In addition, there is a significant contribution to the 
stiffness k from the Neel walls. When the vortex is at 
the edge, the Neel wall energy is clearly lower than when 
it is at the center due to the absence of the parabolic 
segments. Note that there is no change in the energy of 
the straight Neel wall segments, because the total length 
and tension of these segments stay the same as the vortex 
moves. 

The parabolic Neel walls are described by the equa- 
tion {x — X)'^ = {2y ± w){2Y ± w) for the upper and 
lower regions, respectively. By symmetry, the total en- 
ergy of the upper parabolic segment, when the vortex 
position is {X,Y), will be the same as that of the lower 
parabolic segment when the vortex position is {X,—Y). 



Thus, the contribution of each parabolic segment to the 
stiffness constant will be the same. We shall restrict our 
consideration to the lower wall for simplicity. 

It is convenient to derive the energy of the parabolic 
Neel wall by switching to cylindrical coordinates around 
the vortex core. Namely, R — y/{x — X)^ + (y — Y)'^ 
and P = arctan((y — Y)/{x — X)). In these coordinates, 
the lower parabolic wall is given hy R = {w — 2Y)/{1 + 
sin/3). 

The surface tension of a Neel wall is given by Eq. (|C4p . 
where A = \fAf\iQM^ . In the case of the parabolic Neel 
walls, the rotation angle 0o is equal to the angle of the 
wall normal away from the a;-axis. This wall normal 
changes with /3. For the lower wall, (/)o = P/2 + 7r/4. 
The wall hits the edge of the strip when (3 — fioiY) = 
- arcsin((u; + 2Y)/{?,w - 2^)). 

The line element along the parabolic segment is given 
{w - 2y)d/3csc0o/(l -f sin/3). 



by dl = y/dr 



^dl3^ 



Hence, the energy of the lower parabolic Neel wall is 



1-0 

Ep^ = 2V2^ioM'^\t{w-2Y) i 

Jbh 



d(3 



MY) 



1- MP) cot MP) 

l-f sin ^ 

(C14) 

Taking the second derivative of this energy with respect 
to F at F = yields the contribution 



3w 



(C15) 



of the parabolic Neel walls to the spring constant k, where 
the factor of 2 on the left hand side comes from the in- 
clusion of the upper parabolic wall as well as the lower. 

Finally, there is a smaller contribution to the stiffness 
k from the bulk regions of nonuniform magnetization. 
Again, the symmetry of the problem allows us to focus 
on the lower region and simply double our result to find 
the total contribution. Using the cylindrical coordinates 
and definition of /3o, we find that the exchange energy of 
the lower region is given by 

E^^_ = iioM'^t I d^x — 

= ,.M^XH r''^/51n.-(- + 2y) 

J -tt/2 

Jpo 



2ro sin f3 

{w - 2Y) 
ro(l -I- sin/3) ' 



where tq is the radius of the vortex core. If we again 
take the second derivative with respect to Y a.t Y — 0, 
we obtain 



AV2\ 2fio]VPXH 



(C16) 
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In total, we find that 



For t = 20 nm, A = 3.8 nm, and w = 200 nm, this gives 



8A 4A2 



3t 



tw 



2+2 



k = 0.49 



w 



w 



(CIS) 



(C17) 



^ D. A. AUwood, G. Xiong, C. C. Faulkner, D. Atkinson, 

D. Petit, and R. P. Cowburn, Science 309, 1688 (2005). 
^ L. Thomas, M. Hayashi, X. Jiang, R. Moriya, C. Rettner, 

and S. Parkin, Science 315, 1553 (2007). 
^ C. L. Chien, F. Q. Zhu, and J. G. Zhu, Physics Today 60, 

40 (June 2007). 
" G. S. D. Beach, C. Nistor, C. Knutson, M. Tsoi, and J. L. 

Erskine, Nature Mat. 4, 741 (2005). 
^ O. Tchernyshyov and G.-W. Chern, Phys. Rev. Lett. 95, 

197204 (2005). 

® O. A. Tretiakov, D. Clarke, G.-W. Chern, Y. B. Baza- 
liy, and O. Tchernyshyov, Phys. Rev. Lett. 100, 127204 

(2008). 

^ T. L. Gilbert, Phys. Rev. 100, 1243 (1955). 

® T. L. Gilbert, IEEE Trans. Mag. 40, 3443 (2004). 

^ M. J. Donahue and D. G. Porter, Tech. Rep. NISTIR 6376, 

National Institute of Standards and Technology, Gaithers- 

burg, MD (1999), http://math.nist.gov/oommf. 
^° R. D. McMichael and M. J. Donahue, IEEE Trans. Magn. 

33, 4167 (1997). 
" H. Youk, G.-W. Chern, K. Merit, B. Oppenheirner, and 

O. Tchernyshyov, J. Appl. Phys. 99, 08B101 (2006). 

M. R. Freeman, W. K. Hiebert, and A. Stankiewicz, J. 

Appl. Phys. 83, 6217 (1998). 

N. L. Schryer and L. R. Walker, J. Appl. Phys. 45, 5406 
(1974). 

^'^ A. Thiaville and Y. Nakatani, in Spin Dynamics in Con- 
fined Magnetic Structures III (Springer, 2006), vol. 101 of 
Topics in applied physics, pp. 161-205. 
M. Hayashi, L. Thomas, C. Rettner, R. Moriya, and S. S. P. 



Parkin, Nat. Phys. 3, 21 (2007). 

J.-Y. Lee, K.-S. Lee, S. Choi, K. Y. Guslienko, and S.-K. 
Kim, Phys. Rev. B 76, 184408 (2007). 
^'^ K. Y. Guslienko, J.-Y. Lee, and S.-K. Kim (unpublished), 
arXiv:0711.3680vl. 

E. Saitoh, H. Miyajirna, T. Yarnaoka, and G. Tatara, Nar 

ture 432, 203 (2004). 
^® V. W. Doring, Z. Naturforsch 3a, 373 (1948). 
^° G. Tatara and H. Kohno, Phys. Rev. Lett. 92, 086601 

(2004). 

J. Yang, C. Nistor, G. S. D. Beach, and J. L. Erskine, Phys. 
Rev. B. 77, 014413 (2008). 

A. Wachowiak, J. Wiebe, M. Bode, O. Pietzsch, M. Mor- 
genstcrn, and R. Wiesendanger, Science 298, 577 (2002). 

23 K. Y. Guslienko, K.-S. Lee, and S.-K. Kim, Phys. Rev. 
Lett. 100, 027203 (2008). 

2* B. Van Waeyenberge, A. Puzic, H. StoU, K. W. Chou, 
T. Tyliszczak, R. Hertel, M. Fahnle, H. Briickl, K. Rott, 
G. Reiss, et al., Nature (London) 444, 461 (2006). 

2^ S. Komineas, Phys. Rev. Lett. 99, 117202 (2007). 

2® R. Hertel and C. M. Schneider, Phys. Rev. Lett. 97, 177202 
(2006). 

O. A. Tretiakov and O. Tchernyshyov, Phys. Rev. B 75, 
012408 (2007). 

2* A. A. Belavin and A. M. Polyakov, Pis'ma ZheETF 22, 

503 (1975), [JETP Lett. 22, 245 (1975)]. 
2^ G.-W. Chern (unpubhshed). 

3° G.-W. Chern, H. Youk, and O. Tchernyshyov, J. Appl. 
Phys. 99, 08Q505 (2006). 



